Entanglement-based signature of nonlocal dispersion cancellation 
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We derive an inequality bounding the strength of temporal correlations for a pair of light beams 
prepared in a separable state and propagating through dispersive media with opposite signs of 
group velocity dispersion. The presented inequality can be violated by entangled states of light, 
such as photon pairs produced in spontaneous parametric down-conversion. Because the class of 
separable states covers the entire category of classical fields as a particular case, this result provides 
an unambiguously quantum feature of nonlocal dispersion cancellation that cannot be reproduced 
within the classical theory of electromagnetic radiation. 
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I. INTRODUCTION 

The quantum theory of electromagnetic radiation al- 
lows for correlations substantially stronger than those 
permitted in classical models. Such correlations, result- 
ing from the uniquely quantum phenomenon of entan- 
glement, lead to striking effects in different degrees of 
freedom of light including polarization wave vector 
0, and frequency. In the last case, one intriguing ef- 
fect is the preservation of strong temporal correlations 
for fields traveling through dispersive media with oppo- 
site signs of the group velocity dispersion. This effect, 
known as nonlocal dispersion cancellation, has been de- 
scribed theoretically by Franson Q and demonstrated in 
a recent experiment by Baek and coworkers 

The actual role of quantum correlations in nonlocal 
dispersion cancellation has been a subject of a vexatious 
debate. Various classical models have been presented 
to illustrate the gap between predictions of the classi- 
cal and the quantum theories [^-Q- On the other hand, 
analogies between these two theories allow one to repro- 
duce the effect of nonlocal dispersion cancellation using 
classical chaotic light up to certain features, such as ap- 
pearance of a constant background in the detected sig- 
nals In this paper we present an inequality which 
bounds the strength of temporal correlations attainable 
when the fields are prepared in a separable state, i.e. 
when no quantum entanglement between the two light 
beams is present. This inequality, based on a separa- 
bility condition for continuous- variable systems [9|, lldj . 
shows that quantum entanglement is necessary to pre- 
serve all the properties of temporal correlations in non- 
local dispersion cancellation. At the same time, the pre- 
sented result defines limits on how well classical fields can 
reproduce the effect of nonlocal dispersion cancellation. 
This is because the separability condition is satisfied by 
the entire category of classical fields, which correspond 
in the quantum theory to statistical mixtures of coher- 
ent states i.e. are represented by non- negative P- 



ble character. Furthermore, the photodetection theory 
gives exactly the same predictions of photocount proba- 
bilities for classical fields and their quantum mechanical 
counterparts Consequently, the derived inequality 

identifies a quantum feature of nonlocal dispersion can- 
cellation that cannot be mimicked by the classical theory 
of optical radiation. 

This paper is organized as follows. In Sec.|TT]we derive 
the transformation of the propagating fields using the 
Wigner phase space representation. The results are used 
to obtain the inequality in Sec. IIIII Its application to 
a class of time-stationary Gaussian states is analyzed in 
Sec. IIVI Sec. [V] discusses experimental aspects of testing 
the proposed inequality. Finally, Sec. I VII concludes the 
paper. 



II. PROPAGATION 

The physical system under consideration shown in 
Fig. [T] comprises two spatially separated light beams. 
We will describe the fields using annihilation parts of 
frequency-domain electric-field operators £j(j-o) parame- 
terized with the detunings w from the central frequency 
too and satisfying commutation relations 



(1) 



where i,j = 1,2. We will assume that the fields un- 
der consideration have restricted bandwidth so that the 
slowly- varying envelope of the positive-frequency part of 
the electric field operator can be written in the temporal 
domain as: 
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representations permitting only correlations of a separa- 



A transparent way to analyze the evolution of the sys- 
tem is to use the chronocyclic Wigner function [l3j ex- 
tended to a pair of optical beams: 
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WtywMM) = / dT i / dr.e^+^TiigE^ih - %)E^\t 2 - T -f)E { 2 +) (t 2 + f + £)], (3) 



where p is the density matrix characterizing the fields and 
E ( f\t) = [ J E J (+) (t)] t . In the regime of low light intensi- 
ties, the integration of the Wigner function over frequen- 
cies u>i and oj 2 yields the joint probability of detecting 
photocounts at times t\ and t 2 , proportional to the ex- 
pectation value TT[gE{'\t 1 )E^\t 2 )E { 2 +) (t 2 )E{ +) (h)}. 

Propagation through a dispersive medium of length L 
results in multiplying the operators £j(ui) by phase fac- 
tors e lki ^ Uo+u '' L , where kjfuiQ + u) is the wave vector at 
the frequency luq + uj for the jth beam. We will apply 
the standard expansion up to the second order: 

fcj-fwo + u) « fej(wo) + — + ftw 2 , j = 1,2 (4) 

where Uj are group velocities and /3j are parameters char- 
acterizing the group velocity dispersion. A straightfor- 
ward calculation shows that in the Wigner representa- 
tion the effect of dispersive propagation is given by the 
following transformation of the temporal variables: 

t'j =tj + — + 2f5 j Lio J (5) 

V 3 

while the frequencies uij are unaffected. We use prime 
signs to denote the chronocyclic variables after propaga- 
tion. It is worthwhile to note that the above result is 
analogous to the quantum mechanical evolution of a free 
particle [11]. 

III. INEQUALITY 

With Eq. ([5]) at hand, the remaining analysis is ele- 
mentary. Let us consider the time difference between the 
detection events t = t± — t 2 , depicted schematically in 
Fig. Ufa). The statistical properties of this quantity are 
obtained by averaging with the normalized Wigner func- 
tion, which we will denote by angular brackets (...). For 
group velocity dispersion parameters of equal magnitude 
but opposite signs, j3i = —02 = fl, one immediately ob- 
tains that the time difference variance after propagation 
((At') 2 } is given by the following combination of the co- 
variance matrix elements for the initial state of light: 

((Ar') 2 > = ({AT) 2 )+4f3L(ATAn) + {2/3L) 2 ((An) 2 ) (6) 

where f2 = lj\ +uj 2 is the sum frequency for both the pho- 
tons. It is seen that the variance of the time difference is 
modified by two contributions. The first one comes from 
the mixed time- frequency covariance (AtAQ). It can be 
eliminated by an assumption that the quantum statisti- 
cal properties of the beams are invariant with respect to 
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their physical interchange, which flips the sign of r. Al- 
ternatively, suppose that we perform two separate exper- 
iments with swapped positions of the dispersive media, 
and take an arithmetic average ((Ar') 2 ) sym of the two 
variances measured after propagation. As the parameter 
j3 has opposite signs in both the cases, this also removes 
the contribution from the mixed covariance term. Thus 
in the symmetric scenario, the time difference variance is 
simply enhanced by the variance of fi: 

((Ar') 2 ) sym = ((At) 2 } + (2/3L) 2 ((Att) 2 }. (7) 

The question about the role of entanglement in nonlo- 
cal dispersion cancellation now boils down to limitations 
that need to be satisfied by the variances ((At) 2 ) and 
((Ail) 2 ). From the formal point of view, our problem 
is analogous to that of a pair of continuous-variable sys- 
tems, such as particles in one spatial dimension or single 
light modes [1(|. The analogs of r and are, up to 
scaling factors, the difference of positions and the sum 
of momenta of the two particles. It is well known that 
in general the quantum theory allows one to define these 
two observables arbitrarily precisely, which is best illus- 
trated by the celebrated Einstein-Podolsky-Rosen para- 
dox [l5j]. However, if the composite state is separable, i.e. 
it is a statistical mixture of well defined quantum states 
for individual subsystems, then the product of these two 
uncertainties has a lower bound derived by Tan Ex- 
pressed in terms of time and frequency variables, this 
bound takes the form 

((Ar) 2 )((Af7) 2 )>l. (8) 

Inserting this result into Eq. (JTJ) yields: 

((Ar') 2 ) sym >((Ar) 2 ) + g^. (9) 

The above inequality is the central result of this paper. It 
defines the minimum broadening of temporal correlations 
between two light beams during propagation through dis- 
persive media with opposite dispersion signs if no entan- 
glement is present. A violation of this inequality is an 
unambiguous signature that the two beams have been 
initially prepared in an entangled state and it cannot be 
obtained within classical theory of optical radiation. 

IV. GAUSSIAN STATES 

The effect of nonlocal dispersion cancellation has been 
originally discussed for photon pairs generated through 
a decay of pump photons in the process of spontaneous 
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parametric down-conversion [3J. For a monochromatic 
pump, energy conservation means that the frequencies of 
the twin photons must sum up to the frequency of the 
parent pump photon, implying that indeed ((Af2) 2 ) = 0. 
Consequently, temporal correlations are preserved. Such 
photon pairs can be considered as a postselected ensem- 
ble of the complete twin-beam state produced in spon- 
taneous parametric down conversion. In order to expose 
the limitations of the classical theory, it is insightful to 
discuss the Wigner picture of a class of time-stationary 
Gaussian states @, EH that includes the twin-beam state 
as a particular case. This class is characterized by two 
positive spectra of phase- insensitive autocorrelation func- 
tions 

Sj(u) = J dre- iur Tr[gE^\t + T)E^\t)}, j = 1,2 

(10) 

and the complex spectrum of the phase-sensitive cross- 
correlation function 

S#M = / dre™ Tr[gE[ + \t + r)E { 2 + \t)} (11) 



while all the field means and other second-order correla- 
tion functions vanish. The functions Sj(u) can be inter- 
preted as spectral intensity distributions per unit time 
for individual beams (l7j . 



Ti[g£j(ui)£ 3 (u')} = 2ttS ] (u)5{uj - uj'). (12) 



For low intensities, S{ 2 (<*>) yields the probability ampli- 
tude of generating a pair of photons with frequencies ui 
and —uj per unit time, 



TT[g£i{uj)£ 2 {uj')] = 2ttS[1\uj)S(uj+lu'). (13) 



The Wigner function for the Gaussian states under con- 
sideration can be easily found with the help of Wick's 
theorem [l8| to take the form: 



W(ti,wi;t 2> w 2 )«5i(a;i)S2(w 2 )+«(wi+W2) / du e M*i-ta) [S$> (|( Wl - lo 2 + uj)] *S$ (£(wi - uj 2 - vj), (14) 



The first term is simply the product of individual spectra, 
while the second one, confined to the phase space region 
where ui\ + ui 2 — 0, resembles a Wigner function for the 
two-photon probability amplitude s[^(uj) parameterized 
with ti — t 2 and (uii — uj 2 )/2. 

The difference between predictions of the classical and 
the quantum theories lies in the relative magnitude of 
the two terms that form the Wigner function. A simple 
and general way to find relevant relations is to define an 
operator 

J"(A)= [ 2 dw[£i(w) + A£j(-w)] (15) 

J -ei 

where ex, £2 > and A is a complex number, and to 
consider the necessarily nonnegative expectation value 
Tz[qJ-(\)J^ (\j\ > 0. Making use of the commutation 
relations given in Eq. ([1]) and performing frequency inte- 
grals yields 

1 + 5i(w) + 2Re[A*S^V)] + \X\ 2 S 2 (-lj) > (16) 

which holds for an arbitrary complex A. This implies 
that 

\S$(w)\ 2 < [l+5i(w)]S 2 (-w). (17) 

The inequality derived in Eq. (jTT)) shows that in the 
limit of low intensities the marginal spectra Si(oji) and 



S 2 (lo 2 ) are allowed to scale as the square of the two- 
photon probability amplitude. In this case the product 
Si(loi)S 2 (uj 2 ) appearing in the Wigner function calcu- 
lated in Eq. (p~4|) is fourth order in s[ 2 \ui). It can there- 
fore be neglected compared to the second term, which 
exhibits strong time-frequency correlations that are be- 
hind the effect of nonlocal dispersion cancellation and the 
violation of the inequality (J9j) . 

For classical fields, one needs to consider statistical 
averages involving stochastic fields £\ and £ 2 and re- 
peat steps leading to Eq. (117[) . However, the commu- 

tativity of the classical fields implies that IS'^^w)! 2 < 
S\(lS)S 2 {—uj). This means that the product of the 
marginal spectra S\(yi\)S 2 (uj 2 j scales at least quadrati- 
cally with (uj) and contributes to the Wigner function 
with at best the same magnitude as the second term. As 
a result, the temporal distribution of detection events has 
a constant background shown schematically in Fig.[^b), 
which severely affects the time difference variance. For 
time stationary fields it makes ((At) 2 ) diverge to infinity. 
If the source emission time is restricted, for example with 
shutters opened for a finite period, ((At) 2 ) would be en- 
tirely dominated by the constant background. In cither 
case, the inequality ([9]) would be satisfied as expected 
for classical fields. Thus the use of time difference vari- 
ance as a quantitative measure of temporal correlations 
allows one to draw a clear distinction between quantum 
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and classical models of nonlocal dispersion cancellation 
01- 

V. EXPERIMENTAL CONSIDERATIONS 

The results of the experiment conducted by Baek et 
al. [I[ on nonlocal dispersion cancelation are inconclu- 
sive with respect to our criterion, as the reduction in 
the time difference variance has been demonstrated only 
against the case with the negative-dispersion medium ab- 
sent in one arm of the setup. To test the inequality ((9|), 
one needs to measure time difference variance before and 
after dispersive propagation, the latter in a symmetrized 
scenario, and to determine independently effective dis- 
persion. A practical issue in an experimental test may 
be the jitter of single photon detectors that limits the 
timing resolution. In order to assess its impact, let us 
assume that the intrinsic source variance ((Ar) 2 } sourcc 
is enhanced by an additive term ((AT) 2 )j; t ter describ- 
ing the temporal uncertainty of the detector response, 
yielding the actually observed time difference variance 
((Ar) 2 ) obs = {(Ar) 2 ) som . cc + ((Ar) 2 ) jittcl .. For down- 
conversion sources the intrinsic variance ((AT) 2 ) source de- 
pends primarily on the bandwidths of the spectra of in- 
dividual beams, and for broadband down-conversion it 
can be made well below (100 fs) 2 [l9[. In this regime, 
((Ar) 2 )jitter becomes the dominant contribution to the 
observed time difference variance, as the timing resolu- 
tion of currently available single-photon detectors is of 
the order of 50 ps [2p| . 

A prerequisite to observe a signature of entanglement 
is to maintain the violation of the separability criterion 
given in Eq. (J5J even when the temporal correlations are 
affected by the jitter. This gives a condition for the sum 
n of the frequencies of the down-converted photons in the 
form ((Af2) 2 ) < [{(Ar) 2 ) j ittcr]~ 1 , which is effectively a 
constraint on the spectral linewidth of the cw laser pump- 
ing the nonlinear medium. This regime can be reached 
using a narrow linewidth pump laser. Further, in the in- 
equality © both ((At') 2 ) sym on the left-hand side and 
((At) 2 ) are equally affected by the jitter variance, while 
the dispersion-induced term (2/3L) 2 /((At) 2 ) decreases 
when the denominator is enhanced by ((AT) 2 )j;tter- This 



makes it harder to violate an inequality calculated for 
actually observed ((AT') 2 ) sym Q b s and ((Ar) 2 ) b s com- 
pared to that with variances free from jitter effects. 
The remaining issue is the magnitude of the dispersion- 
induced term compared to the time difference variance. 
To make a statistically significant observation, we would 
like the dispersion term to be non-negligible to the ini- 
tial time difference, i.e. (2/?L) 2 /((AT) 2 ) obs « ((AT) 2 ) obs . 
If the observed time difference variance is dominated 
by the detector jitter, this leads to a condition 2/3L ks 
((AT) 2 )jitter- This condition is more challengingto fulfill, 
as in the exemplary experiment of Baek et al. [4| we had 
2/3L rs (8 ps) 2 , which is substantially smaller than typical 
((At) 2 ) jitter- This relation could be improved by trans- 
mitting fields through engineered highly dispersive opti- 
cal fibers [21| and by using detectors based on parametric 
upconversion [22j that enable precise temporal gating by 
combining the signal with a short auxiliary pulse. Over- 
all, a careful choice of the experimental regime should 
enable a violation of the inequality © • 



VI. CONCLUSIONS 

In conclusion, we presented an inequality that reveals 
the role of entanglement in nonlocal dispersion cancella- 
tion. It is worth noting that an analogous discussion can 
be carried out for the spatial degree of freedom when two 
light beams are subjected to suitably arranged diffractive 
propagation (23[. Another interestin g qu estion would be 
to analyze the time interval statistics [24j , whose variance 
could remain finite even in the presence of uniform back- 
ground of detection events that occurs in the classical 
case analyzed above. 
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FIG. 2: Joint time statistics of detection events. Each dot 
(tijtz) represents registration of two photons at times t\ and 
t'2. The object of interest is the marginal distribution of the 
time difference t = ti — t%. In the quantum case (a) all 
the events are temporally correlated, giving a finite value of 
the variance (Ar) 2 . Classically (b), the constant background 
makes the variance divergent. 



FIG. 1: The measurement scheme. Two light beams char- 
acterized by positive- frequency field operators E[ + \ti) and 
E^\t2) propagate through dispersive media of equal lengths 
L but opposite group velocity dispersions /3 and — (3. The 
quantity of interest is the time difference r = ti — ti between 
detection events. 



